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Abstract 

Here we present new adaptive optics observations of the Quaoar-Weywot 
system. With these new observations we determine an improved system 
orbit. Due to a 0.39 day alias that exists in available observations, four pos- 
sible orbital solutions are available with periods of ~ 11.6, ~ 12.0, ~ 12.4, 
and ~ 12.8 days. From the possible orbital solutions, system masses of 
1.3 — 1.5 ± 0.1 x 10 21 kg are found. These observations provide an updated 
density for Quaoar of 2.7 — 5.0 g cm -3 . In all cases, Weywot's orbit is eccen- 
tric, with possible values ~ 0.13 — 0.16. We present a reanalysis of the tidal 
orbital evolution of the Quoaor-Weywot system. We have found that Wey- 
wot has probably evolved to a state of synchronous rotation, and have likely 
preserved their initial inclinations over the age of the Solar system. We find 
that for plausible values of the effective tidal dissipation factor tides produce 
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a very slow evolution of Weywot's eccentricity and semi-major axis. Ac- 
cordingly, it appears that Weywot's eccentricity likely did not tidally evolve 
to its current value from an initially circular orbit. Rather, it seems that 
some other mechanism has raised its eccentricity post-formation, or Weywot 
formed with a non-negligible eccentricity. 
Keywords: Adaptive optics, Kuiper belt, Satellites, dynamics 



1. Introduction 



Large Kuiper Belt objects with diameters D > 1000 km exhi 
range of densities, with values typically larger than 1.5 k cm -3 ( 


}it a broac 
Buie et al.. 


1997: 


Rabinowitz et al.. 


2006; 


Brown et al. 


2010: 


Fraser and Brown. 


20101). 



In addition, many of the larger obje cts are founc 



few percent the size of th e primary (iBrown et al 



2007 



Fraser 



Brown et al 



et al 



2006 



with s mall satellites only a 



2005 



201 Of) and that are 



Brown and Schaller . 



Fraser and Brown 



icy in nature ( iBarkume et al. 



2006 : 



20091 ). A natural explanation is one 



in which during the early phases of planetesimal growth, these large bodies 
accreted a large enough mass sufficiently rapidly to heat up and differen- 
tiate, producing silicate rich cores surrounded by icy mantles. Subsequent 
collisional evolution then stripped predominantly icy material, raising their 
densities above their primordial values. The range of densities exhibited 
by each object would then reflect the relative amounts of collisional bom- 
bardment that object suffered. Another bi-product of this process are the 
satellites of the larger objects which are just collisional fragments from the 
mantle which were not ejected with a high enough velocity to escape, and 
remained bound to their primaries. 
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The large object Qu aoar is an extreme c ase. O bservations of its satellite, 



Weywot, presented by iFraser and Brown! (120101 ) suggest it is the densest 
known KBO. The observed density ofp = 4.2±1.3g cm -3 implies that 



Quaoar may co nsist almost entirely o : 
thin icy veneer (ISchaller and Brownl . 



silic ate material surrounded by a very 



20071 ). Such a result would imply that 



the amount of collisional bombardment suffered by Quaoar was significantly 
higher than that experienced by other large KBOs. 

The Quaoar- Weywot system prese nted a further pecu li arity: Weywot ap- 



peared to be on an eccentric orbit (IFraser and Brownl . |2010| ). This was 
unexpected, as it seemed most likely that tidal evolution would circular- 
ize the orbits of the small satellites on short timescales consistent with 
the satellite of Eris, Dysn omia, which is found on a nearly circular orbit 



( iBrown and Schaller 



20071 ). These two strange properties of the Quaoar- 



Weywot system warrant further investigation. 

Here we present new observations of Quaoar and its satellite Weywot. In 
section [2] we present our observations of this binary system made with the 
Keck 2 telescope and the data reduction steps to identify Weywot within 
the images. In section [3] we present a new determination of Weywot's orbit 
and colour, along with a more accurate determination of the Quaoar- Weywot 
system mass. In section H] we present a re-analysis of the tidal evolution the 
Quaoar- Weywot and Eris-Dysnomia binaries. Specifically, we present eccen- 
tricity and semi-major axis evolution which considers tides raised on both 
bodies and an eccentric orbit of the satellite. In addition, we present some 
tidal evolution simulations of the Quaoar- Weywot system. We discuss pos- 
sible evolution of the orbital parameters of Weywot and Dysnomia. Finally, 
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we finish with a short discussion of the results in section 

2. Observations and Data Reductions 

Observations of the Quaoar-Weywot system were taken on June 7th, 2011 
UT. During that night, Quaoar experienced an appulse with a 10.5" closest 
separation to a R=12.1 magnitude star. This allowed excellent adaptive 
optics image correction in natural guide star mode using the Keck 2 adaptive 
optics facility, resulting in image cores that had Full- Width at Half-Maxima 
of 55-66 milliarcseconds and a Strehl between 0.13 and 0.37, with 10-25% of 
the light in the narrow core. Observations were taken with the Near Infrared 
Camera 2 (NIRC2) in the narrow camera mode resulting in a 0.01" pixel 
scale. Observations were taken in the K' filter using 15 s exposure times and 
a fixed sky position angle was maintained throughout the observations. The 
camera position angle was offset by 0.7 degrees from zero to account for a 
slight rotation that exists between the telescope FOV and the camera. A 3 
point dither pattern was utilized to avoid the bad quadrant of NIRC2 and 
15 images were taken at each dither point. 

Images were reduced with standard techniques. Appropriate darks, bi- 
ases, and dome flats were used to remove instrumental flat-field and bias 
patterns. Median stacks of the 15 images at each dither point were created 
from the de-biased and flat-fielded images. Background levels at a partic- 
ular dither position and in the quadrant containing Quaoar were made by 
averaging the stacks of that quadrant before and after that particular dither. 
The result was a flat image with zero background to within the noise of the 
images. 
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Identification of Weywot was made possible during a portion of the ob- 
servations during which Quaoar passed close to a nearby star. Subtraction of 
Quaoar's point-spread function of sufficient quality for Weywot's easy iden- 
tification was achieved when Quaoar was within 1.2" of the star, roughly 1/4 
of the total of the time in which Quaoar was observed. In addition, the ob- 
servations were taken in fixed sky position angle mode. As a result, the PSF 
and associated speckle pattern rotated by roughly 90° demonstrating that 
the image of Weywot is real and not a PSF artifact. The image subtraction 
results along with a median stack of all images centred on Quaoar's position 
are shown in Figure [TJ Weywot is easily seen South-East of Quaoar. 

To determine accurate astrometric positions of Weywot with respect to 
Quaoar, 4 median stacks with equal equivalent exposure times were produced. 
For each stack, a radial profile of Quaoar's image was subtracted to reveal 
Weywot within the PSF wings. The radial profiles were generated using 15° 
wide radial slices separated by 15° centred on Weywot's position. The resul- 
tant centroids in each image were not significantly altered by large variations 
in the radial slice parameters. The results are presented in Figure [H 

Astrometric positions were determined in three separate ways, PSF match- 
ing using the central 9 pixels of Quaoar's image, using Weywot's photo-centre, 
and gaussian profile fitting, all with a 5 pixel subsampling. The results of all 
three methods were quite consistent, with less than 1 pixel scatter. The final 
accepted positions are the average of all three methods, with a 1 pixel, or 
0.01" uncertainty conservatively adopted to account for possible variations 
due to the centroid and radial profile fitting procedures. Weywot's motion 
can clearly be seen from stack to stack. These measurements are presented 
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in Table [TJ along with available previous astrometric m easurements. Ta 



Fraser and Brown 



blc[T] also includes corrections to the positions presented in 
(120 lOl ) which neglected the cosine of Quaoar's declination at time of the ob- 
servations. 



3. Weywot and its Orbit 

From the four separate stacks, the Quaoar- Weywot system was found to 
have a K'-band flux ratio of 1 : 0.0034 ± 0.0005. Adopting similar albe- 



dos, Weywot has a size ratios o 



observations (IFraser and Brown 



20101) 



in rough agreement with previous 



Weywot's orbit was fit using the same methods as 



Fraser and Brown 



( 120101 ). That is, a maximum likelihood routine was used to consider all 



available detections and non-detections as a result of Weywot falling within 
Quaoar's image. The best-fit orbital solutions are presented in Figure [2] and 
in Table |2j 

The observations do not have a sufficient temporal span to break the 
on-sky mirror degeneracy. In addition, the obs ervations allow an alias i n the 



period of 0.39 days, an alias that was missed by IFraser and Brownl ( 120101 ) who 
only considered the best-fit orbital period of 12.4 days when determining the 
best-fit orbit. Formally, five orbital solutions corresponding to different local 
minima in the likelihood are found including that with the 12.4 day period 
were found when all data are considered. The compatibility of the different 
orbital solutions was tested by comparing the likelihood of each fit to the 
distribution of likelihoods found for that particular orbit. This was done by 
generating random simulated observations from one of the orbital solution. 
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The simulated data were generated with gaussian noise with widths equal to 
the errors in the observed data. The simulated observations were fit and the 
process was repeated to generate a distribution of simulated likelihood values. 
For that fit, the probability P(L ran < L obs ) of finding a random likelihood 
value smaller than the observed value, L obs was determined; values near 1 
indicate a poor solution. From this test, only four periods of nearly ~ 11.68, 
~ 12.04, ~ 12.43, and ~ 12.84 days are possible solutions. Orbits with 
~ 11.3 day periods are excluded at greater than the 99.8% level. The shortest 
acceptable period has P(L ran < L b s ) > 0.8 and the longest has P(L ran < 
Lobs) > 0.95. Thus, while formally consistent with the observations, the 
shortest and longest orbital periods are unlikely the correct ones. Including 
their mirrors, eight acceptable orbits are found. The 12.04 day prograde orbit 
provides the best solution. Additional observations are needed to determine 
which is the true orbit. 

All acceptable orbits include non-zero eccentricities. A circular orbit with 
period P = 12.43 days however, is possible with a moderately acceptable 
Lobs = 8.7. To test the validity of a circular orbit over the best-fit eccen- 
tric orbit, we turn to the likelihood ratio test. We use the ratio to test the 
null hypothesis, that the improvement in the likelihood of an eccentric orbit 
over that of the circular orbit is insignificant. Assuming gaussian distributed 
measurement errors and in the limit of large sample size, the likelihood ratio 
X = — 2 log , Lobs is distributed as a chi-squared distribution with 1 degree of 

L/obs,e— 

freedom. Here, L Q b s is the likelihood of the eccentric orbit fit, and L b s ,e=o is 
the likelihood of the circular restricted fit. Of course, neither of the assump- 
tions about the likelihood ratio necessarily hold true for our data, especially 
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the former; 9 measurements and 2 non-detections certainly do not meet the 
requirement of a large sample size. As a result, we are forced to calibrate the 
likelihood ratio with the use of Monte-Carlo simulations. From the best-fit 
circular orbit we generate an artificial set of observations in the same manner 
as that used to determine P(L ran < L obs ). Both an eccentric and a restricted 
circular orbit were fit to the simulated data, and the likelihood ratio, x was 
determined. This process was repeated to generate a distribution of x val- 
ues. The advantage of this Monte-Carlo process is the complete avoidance 
of the assumption of a large sample size. All of our simulations include the 
same number of simulated observations as we actually have, and the resultant 
range of simulated likelihood ratios will reflect this. 

In our simulations random likelihood ratio values larger than the observed 
value, x = 9-2 were found in less than 1 in 1000 simulations. That is, 
the probability of the null hypothesis, that the eccentric fit results in an 
insignificant improvement over the circular fit, is less than 1 in 1000. This 
demonstrates that the eccentric orbit is a vast improvement over the circular 
fit. Similar improvements were even rarer for other periods and as such, we 
can formally exclude all circular orbit solutions at greater than the 3-a level. 

It should be noted that the uncertainties on the astrometric positions are 
potentially overestimated. This can be seen from the best-fit orbit which 
has P(L ran < L obs ) = 0.14; values near 0.5 are expected. This value suggests 
that the true uncertainty in the astrometric positions of Weywot, on average, 
may be as much as 25% smaller than quoted. It should be made clear that 
this cannot be used to improve the orbital fits as it is virtually impossible to 
know which data points have over estimated uncertainties and which do not. 
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The four possible orbits (and their mirrors) provide a range of masses for 
the Quaoar-Weywot system of 1.3 — 1.4 x 10 21 kg. Given the best estimate 
of Quaoar's diameter of 890 ± 70 km ( IFraser and Brownl . 
density must fall in the range 2.7 — 5.0 g cm -3 . 



20101), Quaoar's 



4. Satellite Orbital Evolution 

The observations confirm that Weywot is on an eccentric orbit with ec- 
centricity e = 0.13 — 0.16. In addition, Weywot is on a highly inclined orbit, 
with an inclination of ~ 15° (152° for the retrograde orbits). The fact that 
Weywot is on an eccentric orbit is surprising given the similar, albeit more 
massive, Eris - pysn omia system which is found on a nearly circular orbit 



( IBrown et al. 



20051 ). 



Here we revisit the question of tidal evolution. In 



( 120101 ) we assumed (as is often done - see for example 



Fraser and Brown 



Noll et al. 



(120081 )) 



that the dissipation from tides raised on the primary is negligible. We have 
found however, that this can be an important effect and here, we consider 
the mutual tides raised on both bodies in determining the orbital evolution 
of the system. 

To get a rough handle on whether or not it is possible that tidal effects 
can explain the differences in the binary orbits between the Quaoar-Weywot 
and Eris-Dysnomia systems, we first consider order of magnitude estimates 
of the tidal evolution timescales in an attempt to determine if reasonable 
tidal solutions exist that allow Weywot to maintain its eccentricity, or have 
it increased to its current value over the age of the Solar syste m. We choose 



to work within the weak friction tidal theory of 



Hut 



(ll98lf ) in which the 
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tidal potential is assume d to 



small time, T. Following 



Hut 



ag the direction of the tide raising body by a 
( )198ll ). we place no restriction on the orbital 
eccentricity, e and expand the equations to first order in inclination, i, which 
has the advantage of decoupling the a— e dynamics from that of i. While little 
is known about the Weywot's inclination with respect to Quaoar's rotation 
pole, we found that the uncertainties in our estimates were largely dominated 
by errors in parameters other than inclination. 

As e and a are the most observationally constrained we consider these 
orbital parameters first. For the mutual tides raised on the binary system, 
the orbit-averaged rate of change of a and e are given by 



< e>- 



and 



■27hr,.{— \ (— 



en 



:i-e 2 ) 



2 xl3/2 



;i+ £>)/3-^(i 



3 2n3/2 



h 
(i) 



< a > 



Qk P T P I — - 1 I — — 



R 



a n 



,2^5/2 



;i + J D)/ 1 -(l-e 2 ) 3/2 /: 



m P j \ a J (1 - e- 
Here, n is the mean-motion of the binary orbit, m, R, T, k with appropriate 
subscripts where P refers to the primary and S refers to the secondary, are the 
masses, radii, tidal time lags, and tidal Love numbers of the two bodies. It is 
useful to rel ate the constant time-lag inherent to the tid al theory utilized here 
19811 ) to t he oft-quoted tidal qua li ty fac tor, Q ( jGoldreich and Peald . 
19661 ). Following lEfroimsky and Laineyl (120071 ) . we express the relationship 
(assuming small lags) as T = (2Q\Q — n\)~ 1 . Whenever quoting Q's, we shall 
use the observed values of flp and n, while taking pseudo-synchronous values 



Hut 
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for Q$ since no observational information is available (note that this choice is 
only sensible if the pseudo-synchronization timescales for the satellites indeed 
turn out to be much shorter than the age of the solar system). The functions 
hi fz-i hi A) an d /s which are used below are functions of eccent ricity whic h 
are of order unity for e C 1 and are presented in closed form by iHutl (ll98ll ) . 
We have defined D as 



D= h^(rn L Y(Rs\ 5 . (3) 
k P Tp \m s J \R P J 1 ; 

Finally, rj = a P mp +™ s (^) 2 (1 - e 2 )~ 1/2 ^f is the ratio of the primary's rota- 
tional angular momentum to the orbital angular momentum where a = jj-^ 
with appropriate subscripts is the specific moment of inertia factor for each 
body. 

Our ultimate goal is to understand why Weywot is currently found on an 
eccentric orbit while Dysnomia is not. Although it is presumed that Weywot 
and Dysnomia were formed through a collisional disruption event, the initial 
orbital and spin configurations of the satellites are poorly constrained. We 
therefore adopt the best-fit parameters to derive order of magnitude estimates 
of the tidal effects at current epoch. In particular, we choose the orbit with 
the highest maximum likelihood (orbit 2 in Tab le 2) for the Quaoar - Weywot 
system, and the orbital parameters found by iBrown et al.l (120051 ) for the 
Eris-Dysnomia system. Because both pos sible rotation per iods of Quaoar 



Ortiz et al. 



(120031 ) greatly exceed 



compatible with the light curve observed by 
the mean motion, here we shall use the shorter period of 8.84 hr and keep 
in mind that the quoted estimates can be translated to the longer period by 



increasing Tp by a factor of 2. For Eris, we adopt a 25.9 hr rotation period 
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as found by iRoe et al.l (120081 ) . 

Relevant uncertain (or unknown) parameters include the radii, rotation 
rates, masses, and moments of inertia, and the effective quality factors of both 
bodies. We estimate radii and masses of the satellites from the observed flux 
ratios of the primary-satellite pair, assuming equal albedos, and take the 
satellite densities to be 1 g cm -3 . The moments of inertia of the bodies are 
assumed to be that o f uniform spheres, while the tidal Love numbers are 



given by ( jPeald . Il999l ) : 



3/2 



k = T^T (4) 



2pgR 

where g is the surface gravity. We adopt for the primary and satellite bodies 
rigidities of fi = 4 x 10 11 and \i — 4 x 10 10 dynes cm -2 appropriate for rocky 
and icy bodies respectively. Note that the variations in density roughly 
cancel the variations in rigidity between the two materials, rendering k an 
approximately unique function of radius. In all cases, adopting a wider range 
of these parameters does not change the main results significantly. Additional 
uncertain parameters will be discussed as needed. Despite the substantial 
number of unknown parameters, we can still determine instructive order of 
magnitude estimates of the tidal evolution of the binary systems at current 
epoch. 

As already mentioned above, the current rotation rates of the satellites 
have not been observed. As will be shown below however, it is possible 
that the satellite's spins have evolved considerably over the age of the solar 
system due to tides. Thus, we adopt pseudo-synchronous angular velocities 
for the satellites to determine estimates of the variation timescales r a = 

a <a> ' 
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r e = — |— . Here, it is important to note that if a satellite is sufficiently 
aspherical, in presence of considerable eccentricity it can become trapped 
in a spin-orbit resonance (jGoldreich and Peald . Il966l ). This phenomenon 
may be of particular importance for Weywot. The associated librations of 
the satellite would provide an additional source of dissipat ion, which is not 



properly taken into account by the tidal theory of lHutl (119811 ). In light of this, 
the degree of dissipation derived here may be viewed as an under-estimate 
for the secondaries. 

The range of r e possible for both systems is presented in Figure [3] as a 
function of the tidal quality factor ratio of the secondary and the primary 
bodies. It should be noted that the range of timescales determined by our or- 
der of magnitude estimates are dominated by the uncertain tidal parameters, 
and the choice of any of the possible orbits for Weywot has little effect. 

From Figure [3] it might be interpreted that the order of magnitude esti- 
mate of the tidal evolution can qualitatively reproduce the contrasting ec- 
centricities of Weywot and Dysnomia. That is, plausible values of ^ can 
be chosen such that the e takes on positive values for Weywot, and negative 
values for Dysnomia. Analysis of Equation [1] reveals that tidal eccentricity 



growth for Weywot requires T S /T P < 40 (corresponding to Qs/Qp 



> I' 

~ 2' 



The eccentricity evolution timescales however, are much too long. Even if 
we set Tg = (which corresponds to a completely non-dissipative secondary, 
and yields the fastest eccentricity growth), the growth timescale for Weywot 
is r c ~ 6.3 x Q P Gyr and the damping timescale for Dysnomia is at least the 
age of the Solar System. These timescales can only be made considerably less 
than the age of the Solar system with unphysically small Q values, certainly 
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incompatible with the values of Q ~ 100 — 500 that have been inferred for 



other icy bo dies such as the Galilean satellites ( IGoldreich and Peale 



Peald . 



1966 



19861 ) . It is important however, to keep in mind that these estimates 



are made using current orbital parameters. 

Using the above argument, we can quantify the sense of tidal evolution 
of the semi-major axis of the Weywot-Quaoar binary. Examination of Equa- 
tion |2] reveals that tidal decay of the orbit requires Ts/Tp > 330 provided 
the current orbit. If we argue that Weywot's eccentricity has been mono- 
tonically increasing due to tidal effects over the system's lifetime, requir- 
ing Ts/Tp < 40, inward tidal migration can be ruled out. The fastest or- 
bital growth however, (again setting T$ = 0) is characterized by a timescale 
r a ~ 14.6 x Q P Gyr, much longer than the age of the Solar system. 

We next turn our attention to the spin-states. The orbit-averaged rate 
of change of the spin velocities, Qp and Qs of primary and secondary bodies 
can be written as 



< tip >= 3k P Tp [ — ] ( — 



m s 



mp 



n 



1 -e 2 ) 6 a P 



e 2 )3/2/5 ^P 

n 



and 



n 



h - (1 - eTVi 



n 



(5) 



(6) 



respectively. The construction of the tidal equations used here is such that 
the spins approach pseudo-synchronization at all times. On its current or- 
bit, Weywot's pseudo-synchronization timescale approximately evaluates to 
tq w ~ 2.4 x 10 -5 x Qs Gyr. Thus, even for the highest quality factors, densi- 
ties, and radii that we can consider, the synchronization timescale for Weywot 
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is no more than 100 Myr. A synchronization timescale roughly a factor of 
50% shorter is found for Dysnomia. As a result, it seems likely that both 
Weywot and Dysnomia will have attained states of near synchronous rota- 
tion by the current epoch, although we again stress that in reality, Weywot's 
significant eccentricity may allow for capture into a high-order spin-orbit 
resonance. 

Comparison of Equations [5] and [6] reveals that the synchronization timescales 
for the primaries is a factor of ~ (jjj^) (^j 3 (g) 3 ~ > 10 5 longer 

than for the secondaries. Thus, the synchronization timescales for the primary 
bodies are > 100 Gyr. As a consequence, the observed rotation rates of both 
Quaoar and Eris are likely to be very close to their primordial values. 

Finally, we consider the binary orbital inclination. That is, the incli- 
nation, i, of the orbit with respect to the primary's spin axis. For small 
inclinations, the variation is well approximated by 



< i >= -3k P Tp I — ] ( — 



mg \ / Rp \ 5 n 2 i 



m P ) V a ) (1 - e 2 ) 13/2 



f 2 1/1 \ fip / 2 \3/2 



-_± i-ip(l- e ^/a 



(7) 

For Weywot near a state of synchronous rotation, the inclination damping 
timescale is Tj Q = 70 x Q P Gyr for Quaoar. A similar timescale is 

found for Eris. Even if the tidal damping factors are extraordinarily small, 
the inclination damping timescale still remains longer than ~ 700 Gyr. Thus, 
while little constraint is available as to the actual values of % for either system, 
it seems likely that both systems have preserved their primordial inclinations 
to the current day. It is noteworthy that the spin-axes of the secondaries align 
orders of magnitude faster than their primaries. 

15 



Our findings suggest that at present, the tidal evolution of the Weywot- 
Quaoar binary is generally slow. This may not have always been the case. 
Consequently, it is important to consider evolutionary paths of the system, 
with an eye towards determining the conditions under which Weywot can 
attain its eccentricity over a period of 4.5 Gyr. To address this, we consider 
tidal evolution simulations with the current observed state as a required end 
result of the simulations. Simulations were run with particular initial condi- 
tions searched to produce the required end state after 4.5 Gyr of integration. 
For these simulations, again we work within the context of the constant time- 
lag theory presented above and envision a giant impact formation scenario 
of the binary, similar to that inferred for the Moon. That is, we assume 
Weywot to have emerged from a proto-satellite disk on a near-circular orbit 
which then migrates outward achieving its current eccentricity 4.5 Gyr later. 

Taking Quaoar's rotation period to be 8.84 hrs, such an evolutionary se- 
quence could not be found with reasonable values of Q. More precisely, for 
an initial eccentricity of eo = 0.01 and a dissipation factor Qw = 100, an 
evolution that produced the observed eccentricity requires unphysically small 
values of the dissipation factor for Quaoar, Qq < 1. Recall that physically, 
Q = 1 implies complete dissipation of all the energy stored in a single libra- 
tion cycle (it should also be noted that the weak-friction theory utilized in 
this work is strictly speaking inapplicable when Q < 10.) Adjusting Qw to 
different values did not improve our results. No satisfactory result could be 
found with physically reasonable values of Qq. 

The reader should be reminded that in our integrations, we assumed 
Weywot had an initial inclination sufficiently low such that cos? ~ 1. This 
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assumption may be invalid; a non-zero initial inclination might be possible. 
The effect of a moderate inclination on the eccentricity evolution arises only 
indirectly through the change in semi-major axis (at second order). As a 
result, even inclinations cannot change our results. It seems that simple tidal 
evolution is insufficient to drive Weywot's eccentricity to its current value over 
the age of the Solar system and it appears that tidal eccentricity excitation is 
not the dominant mechanism responsible for Weywot's orbit. The origins of 
this discrepancy could potentially be attributed to the simplicity of the tidal 
model we have utilized. In particular, we have ignored the possibly complex 
frequency and amplitude dependence of the tidal quality factors, a possibility 
that is difficult to test without a w ell-grounded theory of tidal dissipation 

for an excellent review). 



sec 



Efroimsky and Williams 



2009 



A plausible alternative scenario is one where Weywot's eccentricity is ex- 
cited by some other mechanism and is then maintained due to a very long 
tidal circularization timescale. After all, as shown above, maintenance or 
even slow growth of eccentricity is naturally attained at the current orbit, 
given reasonable tidal quality factors. Possible mechanisms include collisions 
and close passages with other large bodies, or resonance passage during oth- 
erwise smooth tidal evolution. This would imply the presence of additional 
bodies in the system. It is further possible that Weywot did not coalesce 
out of a debris disk in orbit about Quaoar, but rather formed as a coherent 
collisional fragment that was ejected onto an orbit with non-zero eccentricity. 
The nature of such possibilities however, remains elusive. 
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5. Discussion 

The observations presented here have significantly improved the accu- 
racy of Weywot's orbital parameters given a particular orbital period and 
orientation. All available data however, are still unable to break the mir- 
ror degeneracy, or remove the aliasing which allows orbital solutions with 
four different periods. Regardless of the correct orbit, the observations have 
confirmed Weywot's significant orbital eccentricity, requiring an eccentricity 
larger than 0.1 to explain the observations. 

Analysis of the Quaoar-Weywot tidal evolution reveals that tidal dissipa- 
tion in Quaoar is not necessarily negligible. It appears however, that tidal 
evolution is not the dominant mechanism responsible for Weywot's eccentric- 
ity. A plausible history is one in which Weywot primordially low eccentricity 
was driven to one by some other mechanism, and subsequent tidal evolution 
had only small part in producing its current orbit. Such mechanisms might 
include collisions, or near passages with other massive bodies. Alternatively, 
it is possible that Weywot formed with a non-zero eccentricity, hinting that 
Weywot might be a coherent collisional fragment rather than a body which 
coalesced out of a disk of material in orbit about Quaoar. 

We contrast this with the Eris-Dysnomia system which exhibits a simi- 
larly slow tidal evolution. This suggests that what ever mechanism produced 
Weywot's eccentricity eg. formation with high-e, was avoided by Dysno- 
mia. This presents an intriguing diversity in small satellites around the large 
Kuiper Belt Objects. 

Our observations have confirmed Quaoar's large mass, which must be 
larger than 1.3 ± 0.1 x 10 21 kg. Despite the remaining uncertainty in its 
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mass, and the large uncertainty in its diameter, Quaoar still appears to be 
the densest known Kuiper belt object with a range of densities of p = 2.7 — 
5.0 g cm -3 . 

The lower limit on Quaoar's density places it among the class of large, 
collisionally disrupted KBOs, along with objects like Haumea. These objects 
have apparently suffered significant post-differentiation, collisional evolution 
which has stripped away the majority of these object's icy content leaving 
only a thin veneer of ice left on their surfaces. In the case of Haumea, a 
single nearly catastrophic co llision appears to be the culprit (see for ex- 



ample 



Leinhardt et al 



20101 ). Quaoar's unusually high density suggests 



that it may be the result of a different collisional scenario than that which 
disrupted Haumea. On one hand, if a single collision is responsible, then 
Quaoar may have suffered a much more head-on impact than did Haumea, 
resulting in a much higher loss of its icy mantle. This may also explain 



(Ortiz et al. 


2003 


) - as compared 


g Haumea (] 


-labinowitz et al.. 


2006 



Lacerda et al. 



ongated and 



20081 ). Alter- 



natively, Quaoar may have suffered a so-called hit-and-run impact in which 
Quaoar collided with a m uch larger body dest roying Quaoar's icy mantle, 
but leaving its core intact ( lAsphaug et all 120061 ). On the other hand, such a 
high density may be the result of multiple impacts which have slowly chipped 
away at Quaoar's primordial icy mantle. 
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Table 1: Quaoar and Satellite Positions. A RA and 
A DEC are the differences in position between the 
Quaoar and the s atellit e. The first 7 rows are from 
Fraser and Brown ( 201ol ). 



Epoch Satellite-Quaoar Offsets 1 



(JD+2453000) 


A R.A. (arcsec.) 


A Dec. (arcsec.) 


781.38031 


0.328 ±0.01 


-0.119 ±0.01 


1179.12990 


0.303 ±0.03 


-0.135 ±0.03 


1535.70263 


-0.49 ± 0.04 


-0.02 ±0.04 


1540.56061 


0.34 ±0.04 


-0.08 ±0.04 


1546.18353 


-0.45 ± 0.04 


0.09 ±0.04 


1550.31485 






1556.44075 






2719.82856 


0.17 ±0.01 


-0.16 ±0.01 


2719.89737 


0.18 ±0.01 


-0.16 ±0.01 


2719.95760 


0.19 ±0.01 


-0.17 ±0.01 


2720.01759 


0.20 ±0.01 


-0.16 ±0.01 



1 Values quoted in lFraser and Brownl (|2010f ) corrected to include cos Dec term 
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Table 2: Orbital Solutions 1 





Orbit 1 


Orbit 2 


Orbit 3 


Orbit 4 


log(-L b s ) 


7.40 


3.55 


4.1 


10.2 




0.74 


0.14 


0.24 


0.04 


Period (days) 


11.6836 ±0.0002 


12.0457 ± 0.0002 


12.4314 ± 0.0002 


12.8428 ± 0.0002 


Semi-major Axis (10 4 km) 


1.31 ±0.02 


1.35 ±0.02 


1.39 ±0.02 


1.45 ±0.02 


Eccentricity 


0.161 ±0.006 


0.152 ±0.005 


0.137 ±0.006 


0.164 ±0.009 


Inclination (dcg) 


14.8 ±0.7 


14.8 ±0.7 


15.8 ±0.7 


15.5 ±0.7 


Longitude of Ascending Node (deg) 


15.5 ±0.8 


2.43 ±0.7 


1.0 ±0.7 


1.7 ±0.8 


Argument of Perihelion (deg) 


358.6 ±0.8 


8.44 ±0.7 


335.0 ±0.7 


301.8 ± 0.7 


Epoch of Perihelion -2454551 (JD) 


2.67 ±0.1 


2.64 ±0.1 


1.58 ±0.1 


0.658 ± 0.1 


System Mass (10 21 kg) 


1.30 ±0.09 


1.34 ±0.09 


1.39 ±0.09 


1.47 ±0.1 


Density (g cm -3 ) 2 


3.5 ± 0.8 


3.6 ±0.9 


3.8 ±0.9 


4.0 ±0.9 


Mirrors 


log(-i D 6 s ) 


8.71 


3.53 


3.77 


10.2 




0.83 


0.13 


0.18 


0.04 


Period (days) 


11.6856 ±0.02 


12.0476 ± 0.0002 


12.4331 ± 0.0002 


12.8443 ± 0.0002 


Semi-major Axis (10 4 km) 


1.30 ±0.02 


1.37 ±0.02 


1.41 ±0.02 


1.46 ± 0.02 


Eccentricity 


0.154 ±0.006 


0.146 ±0.006 


0.128 ±0.005 


0.157 ±0.008 


Inclination (dcg) 


151.2 ±0.7 


152.1 ±0.7 


150.7 ±0.7 


150.5 ± 0.7 


Longitude of Ascending Node (deg) 


330.7 ±0.9 


341.0 ±0.8 


341.3 ±0.8 


339.2 ± 0.8 


Argument of Perihelion (deg) 


3.5 ± 0.8 


22.5 ±0.7 


344.1 ±0.7 


311.3 ± 0.7 


Epoch of Perihelion -2454551 (JD) 


2.5 ± 0.1 


2.9 ±0.1 


2.7± 0.1 


0.7 ±0.1 


System Mass (10 21 kg) 


1.28 ±0.09 


1.41 ±0.09 


1.43 ±0.1 


1.51 ±0.1 


Density (g cm -3 ) 2 


3.5 ± 0.9 


3.8 ±0.9 


3.9 ±0.9 


4.1 ±0.9 



1 Uncertainties are the 1 — a confidence limits on the relevant parameter with all others 

held at their best-fit values. Angles are with respect to the J2000 ecliptic. 
2 Assuming a diameter for Quaoar of D — 890 ± 70 km. 
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Figure 1: a: Median stack of all images centred on Quaoar. b: Median stack of first 
quarter of the images. The radial profile of Quaoar's image has been removed in the lower 
quadrant containing Weywot. c: As in b, but for the second quarter of the images, d: 
Stack of all images for which PSF subtraction was possible, e: As in b, but for the third 
quarter of the images, f: As in b, but for the fourth quarter of the images. In all images, 
North is up, East is to the left. 
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Figure 2: Weywot's observed offsets from Quaoar (shown at figure centre) and best- 
fitting orbit. Crosses are the observed position along with uncertainty at each epoch 
(see Table [IJ. The best- fitting orbit (orbit 2 from Table [2]) is shown as an ellipse, and 
the predicted positions at each observation from that orbit are shown as small circles. 
The large thick-lined circle represents the region in which Weywot would not have been 
detected in previous observations (see Fraser and Brown, 2010 for details). 
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Figure 3: Eccentricity evolution timescales versus the ratio of tidal quality factor of the 
secondary and the primary bodies. The timescales are shown with Qp = 100 and scale 
linearly with Qp . Black solid and dashed lines correspond to the Quaoar-Weywot s ystem 



Ortiz et al 



((20030 . Red 



with Quaoar rotating with 8.64 and 17.68 hr periods presented by 1 
dash-dotted line shows the timescales for the Eris-Dysnomia system. Here we have adopted 
densities for the satellite of 1 g cm -3 and a radius and albedo for Eris of R p = 1170 km 
and 0.96 (Sicardy, B. personal communication). 
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